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Abstract 

The decay of overlap between a wave packet evolved with a Hamiltonian Ti and the 
same state evolved with TC+T, serves as a measure of the decoherence timer^. Recent 
experimental and analytical evidence on classically chaotic systems suggest that, 
under certain conditions, depends on TL but not on E. By solving numerically a 
Hamiltonian model we find evidence of that property provided that the system shows 
a Wigner-Dyson spectrum (which defines quantum chaos) and the perturbation 
exceeds a crytical value defined by the parametric correlations of the spectra. 



The existence of chaos in classical mechanics is manifested in the evolution of 
a state as an extreme sensitivity to the initial conditions. Quantum mechanics, 
on the opposite, does not show this sensitivity.. This has raised several prob- 
lems in a dynamical definition of quantum chaos. In particular, numerical [1] 
and experimental [2] studies show that time reversal can be achieved with great 
accuracy. Therefore, the search for a quantum definition of chaos, lead to inves- 
tigate the spectral properties [3] of quantum systems whose classical equivalent 
is chaotic. Quantum chaos appears as the regime in which the properties of 
the eigenstates follow the predictions of the Random Matrix Theory (RMT). 
In particular the normalized spacing between energy levels s = (£i+i — £i)/Ae, 
with Ae the mean level spacing, should have a probability distribution given 
by the Wigner Dyson distribution P\yd( s ) = ( 7TS /^) exp(— tcs 2 /A) for an or- 
thogonal ensemble and Pwd( s ) = (32s 2 /7r 2 ) exp(— 4s 2 /n) for the unitary 
ensemble. 

An infinite set of interacting spins is an example of a many-body system 
which is chaotic in its classical version (lattice gas) and hence it is expected 
to present the quantum signatures of chaos in the spectrum. The dynamics of 
this particular system can be studied by Nuclear Magnetic Resonance (NMR). 
Surprisingly, the "diffusive" dynamics of a local excitation exp[— i7-it R ]|0 > can 
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be reversed [4], generating a Polarization Echo M at time 2t#. In this case 
Ji is the many-body Hamiltonian of a network of spins with dipolar interac- 
tion. To accomplish this, the transformation Ti — > — (7i + £) at time t R is 
performed with standard NMR techniques [4] . This transformation is possible 
due to the anisotropic nature of the dipolar interaction. The perturbation E 
is a non-invertible component of the Hamiltonian. In some systems, the only 
contribution to E is proportional to the inverse of the radio frequency power 
and hence can be made arbitrarily small. In one body systems the Polariza- 
tion Echo (i.e. magnitude of the excitation recovered at time 2t R ) can then be 
written exactly as: 

M(t) = | (0| exp[i(H + Y)t/?i} exp[-iHt/h] |0) | 2 , (1) 

where |0) is the initial wave function, Ti the unperturbed Hamiltonian and 
E the perturbation which can be associated with an environmental distur- 
bance. Then the magnitude of experimental interest is the overlap between 
the same initial wave function evolved with the two different Hamiltonians, 
H and — {TL + E). We should note that the second evolution can be seen as 
an " imperfect time reversal" of the wave function!. Consistently, more than 
10 years ago Peres [5] proposed that dynamical signatures of quantum chaos 
should be searched on the sensitivity to perturbations in the Hamiltonian. 
Actually, for a classically chaotic system a perturbation in the initial con- 
ditions is equivalent to a perturbation in the Hamiltonian. The experiments 
show that M decays rapidly with tn with Gaussian law [6] indicating a pro- 
gressive failure in rebuilding the original state. We can define a decoherence 
time T0 from this failure, as the width of the Gaussian. This is found[7] to be 
roughly independent of £ and it extrapolates to a finite value when E — > 0. 
Using a semiclassical one-body analytical approach in classically chaotic sys- 
tems characterized by a Lyapunov exponent A, Pastawski and Jalabert[8] have 
shown that there is a regime where the attenuation of M it is independent of 
the perturbation E and becoming = A. This non-perturbative result is 
valid for long times and as long as E does not change the Hamiltonian nature. 
Our general aim is to find numerical evidence of this regime where M(t) is 
independent of E considering the simplest Hamiltonian that could model spin 
diffusion. 

In this work, we study one-body Hamiltonian in quasi-lD systems with N 
states which we called the Star's necklace model. More specifically, we use a 
tight binding model of a ring-shaped lattice with on-site disorder, hopping 
matrix elements V, and a magnetic flux <3> perpendicular to the plane of the 
ring (see inset in Fig.l). Let us discuss the general features through one rep- 
resentative class, each star has 20 sites and there are L = 35 beads in the 
necklace which makes N = 700. In our case, perturbation acts only between 
two star beads: E(5$) = V exp[i27r$/$ ](exp[i27r<5$/$ ] - 1) |1) (L\+c.c. Bra 
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and ket states contain orbitals in the star. The localized wave packet with 
energy (0| Ti |0) ~ moves along the string and contains only 3N/4 states. 
Anderson disorder is W = 3V which gives a r imp ~ 3h/V and $ = 0.1$ - We 
verify that for <5$ = the dynamics of the system follows a diffusive law, and 
that for all 5$ the statistics of the eigenvalues correspond to those predicted 
by RMT (see Fig.l). We interpret these fact as a signature of chaos. How- 
ever, when studying the parametric correlations of the energy spectrum[9] as 
a function of 5$ the correlation function is definite positive. A critical value 
5$ c ~ 0.1 $ c a n be extrapolated from the strong decay. 

For small t the decay of M(t) is Gaussian like with a characteristic time scaling 
with E exactly as could be predicted by perturbation theory[5]. After a certain 
time of the order of the collision time r imp it becomes a stretched exponential 
with a characteristic time independent of E, 



with v ~ 0.87. and ~ 18h/V (see Fig. 2). The constant arises from 
finite size effects. Nonetheless, if the perturbation does not exceed the critical 
threshold consistent with that of the correlation function, increases with 
decreasing E. On the other hand, if the perturbation is large <5$ > 5$ c , one 
gets Moo ~ l/N. (see inset of Fig. 2). The fact that E c goes to zero when N 
goes to infinity [9], together with results of Fig. 2 for a finite system, could be a 
signature of the existence of a nontrivial thermodynamic limit lim lim ¥«, = 



results on the variation of with the amount of disorder and system size are 
consistent with this hypothesis. 

In order to present in a graphical way the physical phenomena of decoherence, 
we calculated the weight of the wave functions evolved with the two different 
Hamiltonians and the overlap between them as a function of the layer. The re- 
sults are shown in Fig. 3. It can be seen that the evolution of the probabilities 
described by the perturbed and unperturbed Hamiltonians are not signifi- 
cantly different, i.e. they would give the same coarse grained values. However, 
the perturbation causes spatial fluctuations in the phases which produce an 
integral overlap that decays to zero. 

To sum up, our numerical calculations of the polarization echo M in a simple 
chaotic model indicate that there is a regime of the perturbation where the 
decoherence time depends only on the perturbation. This basic feature is also 
found in experiments [7] and in other theoretical models. Some differences in 
the details remain, such as the value in the exponent v. According to prelim- 
inary numerical evidence v could be related to the particular topology [10] in- 
duced by the matrix elements in the Hamiltonian model. Different systems 



M(t) ~ expH/r^ + M c 



(2) 




E-.0JV^oo 

= 1. Preliminary 
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Fig. 1. Probability distribution of the normalized energy level spacing s for the 
system studied. The black line is the Wigner-Dyson distribution for the unitary 
ensemble, which is in reasonable agreement with the numerical data. In the inset is 
shown an schematics of the Hamiltonian model, where the system has the shape of 
a ring, and the in-layer sites are fully connected. 

Fig. 2. Plot of the decoherence time of the system according to 2 as a function of 
<5<E>. The curve is for a system with 20 sites per layer and a diameter of 35 sites, the 
amount of disorder is W = 3. The inset shows the dependence of the asymptotic 
value Mqo with M>. The straight line corresponds to 4/3N. As it can be seen 
reaches this value when 5<& > (5$ c 

such as disordered cylinders, maximally connected Hamiltonians (RMT) and 
chaotic stadiums should be studied in order to characterize this dependence. 
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Fig. 3. This figure shows that while the wave functions evolved with 7i and Tt + £ 
have very similar probability densities, the overlap between them decays rapidly to 
zero due to interference effects. In thick solid line, profile of the initial wave function 
and overlap of the two initial wave functions as a function of the layer respectively. In 
thin solid lines, profile of the wave functions evolved with two different Hamiltonians 
as a function of the layer. Thick dotted line, overlap between the two evolved wave 
functions. The system has 15 sites per layer and a diameter of 35 sites, W = 3, 
t = 2h/V, and 5& = 0.1 <&q. The sum over layers for this overlap is equal to 1.4 
10~ 2 
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S = 8e / Ae 
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